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Abstract 

We compute the imaginary part of the 2-loop vertex corrections in the QCD Fac- 
torization framework for hadronic two-body decays as B ^ vrvr. This completes 
the NNLO calculation of the imaginary part of the topological tree amplitudes 
and represents an important step towards a NNLO prediction of direct CP asym- 
metries in QCD Factorization. Concerning the technical aspects, we find that 
soft and collinear infrared divergences cancel in the hard-scattering kernels which 
demonstrates factorization at the 2-loop order. All results are obtained analyti- 
cally including the dependence on the charm quark mass. The numerical impact 
of the NNLO corrections is found to be significant, in particular they lead to an 
enhancement of the strong phase of the colour-suppressed tree amplitude. 
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1 Introduction 



Charmless hadronic B decays provide important information on the unitarity triangle 
which may help to reveal the nature of flavour mixing and CP violation. In order to 
exploit the rich amount of data that is currently being collected at the B factories, a 
quantitative control of the underlying strong-interaction effects is highly desirable. In the 
QCD Factorization framework [T] the hadronic matrix elements of the operators in the 
effective weak Hamiltonian simplify considerably in the heavy-quark limit. Schematically, 



where the non-perturbative strong-interaction effects are encoded in a form factor ^ 
at = 0, decay constants /m and light-cone distribution amphtudes (pM- The short- 
distance kernels T/ = 0{1) and T/^ = 0{as) provide the basis for a systematic imple- 
mentation of radiative corrections; the former contain the short- distance interactions that 
do not involve the spectator antiquark from the decaying B meson [vertex corrections) 
and the latter describe the ones with the spectator antiquark {spectator scattering) . 

The next-to-leading order (NLO) corrections to the kernels T/'^^, which constitute an 
0{as) correction to naive factorization, are already known from [T]. Recently, the next- 
to- next-to-leading order (NNLO) corrections to T/^ have been computed [21 El IH El E] • 
Due to the interaction with the soft spectator antiquark, the spectator scattering term 
receives contributions from the hard scale ~ mi, and from an intermediate (hard-coUinear) 
scale ~ (AqcD^^fe)^''^- Both types of contributions are now available at (1-loop), 
indicating that the NNLO corrections are numerically important. 

In this work we compute NNLO corrections to T/ for the so-called topological tree 
amplitudes (which arise from the insertion of current- current operators). In contrast to 
the spectator scattering term, the vertex corrections are dominated solely by hard effects 
and amount to a 2-loop calculation. In particular, we address the imaginary part of the 
hard-scattering kernels which is the origin of a strong rescattering phase shift that blurs 
the information on the weak phases. As an imaginary part is first generated at 0{as), 
higher order perturbative corrections are expected to significantly influence the pattern of 
strong phases and hence direct CP asymmetries. Our calculation represents an important 
step towards a NNLO prediction of direct CP asymmetries in QCD Factorization. 

The outline of this paper is as follows: In Section [2] we present our strategy for 
the calculation of the topological tree amplitudes by introducing two different operator 
bases. Section El contains the technical aspects of the 2-loop calculation. In Section H] 
we show how to extract the hard-scattering kernels from the hadronic matrix elements. 
Our analytical results can be found in Section El The numerical impact of the NNLO 
vertex corrections is investigated in Section El and we finally conclude in Section [71 A 
more detailed presentation of the considered calculation can be found in [7]. 
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2 Choice of operator basis 



In view of the calculation of topological tree amplitudes, we restrict our attention to the 
current-current operators of the effective weak Hamiltonian for b —>■ u transitions 

Kff = ^ VubV:, (CiQi + C2Q2) + h.c. (2) 

Due to the fact that we work within Dimensional Regularization (DR), we also have to 
consider evanescent operators [H] • These non-physical operators vanish in four dimensions 
but contribute at intermediate steps of the calculation in d = 4 — 2e dimensions. As the 
imaginary part has effectively 1-loop complexity with respect to renormalization at (9(q;^), 
the considered calculation only requires 1-loop evanescent operators. For our purposes 
the complete operator basis is thus given by 

Qi = [uiYLbi] [dj'j^Luj] , 
Q2 = [ua^L bj] [dj-ff,L Ui] , 

El = [uiYYYLbi] [dj-f^^,-fpLuj] - (16 - 4e) Qi, 

E2 = [u^i^l^YL bj] [d^i^^uipL u,] - (16 - As) Q2, (3) 

where i,j are colour indices and L = 1 — 75. The operator basis in ([3]) has been used 
in all previous calculations within QCD Factorization [U, El |3l H] . We refer to this basis 
as the traditional basis for convenience and denote the corresponding Wilson coefficients 
and operators with a tilde. 

It has been argued by Chetyrkin, Misiak and Miinz (CMM) that one should use 
a different operator basis in order to perform multi-loop calculations [S]. Although the 
deeper reason is related to the penguin operators which we do not consider here, we prefer 
to introduce the CMM basis in view of future extensions of our work. This basis allows to 
consistently use DR with a naive anticommuting 75 to all orders in perturbation theory. 
In the CMM basis the current-current operators and corresponding 1-loop evanescent 
operators read (indicated by a hat) 

Qi= [u,rLT^pj] [d,j,LT^iUi], 
Q2 = [uiY'Lbi] [dj'y^Luj] , 

E, = [ua'-YYLT^pj] [dkl.l.lpLT^^ui] - 16 Qi, 

E2 = [ua^YYLbi] [dj^^^^^pLuj] - 16 Q2, (4) 

with colour matrices and colour indices i,j, k, I. 

Comparing the operator bases in ([3]) and (jl]) we observe two differences: First, the 
two bases use different colour decompositions which is a rather trivial point. More im- 
portantly, they contain slightly different definitions of evanescent operators. Whereas the 
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Figure 1: Generic 1-loop diagram with different insertions of a four- quark operator 
Qi . The upper lines go into the emitted meson M2 , the quark to the right of the vertex 
and the spectator antiquark in the B meson (not drawn) form the recoil meson Mi. 

definitions in the CMM basis correspond to the simplest prescription to define evanescent 
operators, sub leading terms of 0{e) appear in the one of the traditional basis. These 
terms have been properly adjusted such that Fierz symmetry holds to 1-loop order in d 
dimensions. 

We follow the notation of [ID] and express the hadronic matrix elements of the effective 
weak Hamiltonian in terms of topological amplitudes aj(MiM2). E.g. the tx^tx^ 
decay amplitude is written as 



where At,t, = iGp/V^ "^s-^+'^(0)/7r. The amplitude ai(MiM2) is called the colour- 
allowed tree amplitude which corresponds to the flavour content [qgh] of the decaying 
B meson, [qgu] of the recoil meson Mi and [u(l\ of the emitted meson M2. The colour- 
suppressed tree amplitude a2{MiM2) then belongs to the flavour contents [g^fe], [qsd] and 
[mm], respectively. For more details concerning the deflnition of the topological amplitudes 
we refer to section 2.2 in [TO] . 

According to this deflnition, the left (right) diagram in Figure [T] contributes to the 
tree amplitude ai (02) • On the technical level these two insertions of a four-quark 
operator correspond to two different calculations. Instead of performing both calculations 
explicitly, we may alternatively compute the amplitude 02 by inserting Fierz reordered 
operators into the left diagram of Figure [H To do so, it is essential to work with an 
operator basis that respects Fierz symmetry in d dimensions. As we have argued above, 
Fierz symmetry indeed holds to 1-loop order in the traditional basis from ([3]) which allows 
us to derive 02 directly from ai. 

We conclude that the CMM basis is the appropriate choice for a 2-loop calcula- 
tion whereas the traditional basis provides a short-cut for the derivation of the colour- 
suppressed amplitude. We therefore pursue the following strategy for the calculation of 
the imaginary part of the NNLO vertex corrections: We perform the explicit 2-loop cal- 
culation in the CMM basis using the flrst type of operator insertion in the left diagram 
from Figure [H From this we obtain ai{Ci). We then transform this expression into 
the traditional basis which yields ai(Ci) and flnally apply Fierz symmetry arguments to 
derive 0^2 (Q) from ai(C*j) by simply exchanging Ci ^ C2. 
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3 2-loop calculation 



The core of the considered calculation consists in the computation of the matrix elements 



to 0{a1) which represents a 2-loop calculation. As will be described in Section W2[ only 
(naively) non-factorizable diagrams with at least one gluon connecting the two currents in 
the left diagram of Figure [T] have to be considered here. The full NNLO calculation thus 
involves the 2-loop diagrams depicted in Figure [21 but only about half of the diagrams 
give rise to an imaginary part. It is an easy task to identify this subset of diagrams since 
the generation of an imaginary part is always related to final state interactions. 

In our calculation we treat the partons on-shell and write qi = uq, q2 = uq and 
p' = p — q satisfying = and p^ = 2p ■ q = ml (with u = 1 — u). We use DR for the 
regularization of ultraviolet (UV) and infrared (IR) divergences and an anticommuting 
75 according to the NDR scheme. We stress that we do not perform any projection onto 
the bound states in the partonic calculation. We instead treat the two currents in the 
four-quark operators independently and make use of the equations of motion in order to 
simplify the Dirac structures of the diagrams. In order to calculate the large number of 
2-loop integrals we proceed as follows: Using a general tensor decomposition of the loop 
integrals, we essentially deal with the calculation of scalar integrals. With the help of 
an automatized reduction algorithm, we are able to express several thousands of scalar 
integrals in terms of a small set of so-called Master Integrals (Mis). The most difficult 
part finally consists in the calculation of these Mis. 

In the remainder of this section we present some techniques that we have found useful 
for the considered calculation. We sketch the basic ideas of the aforementioned reduction 
algorithm and discuss several techniques for the calculation of the Mis. We refer to the 
references quoted in the following subsections for more detailed descriptions (see also [7]). 

3.1 Reduction to Master Integrals 

Any scalar 2-loop integral in our calculation can be expressed as 



where the Si are scalar products of a loop momentum with an external momentum or of 
two loop momenta. The Vi denote the denominators of propagators and the exponents 
fulfil ni,mi > 0. The scalar integrals themselves depend on the convolution variable u 
in the factorization formula ([T]). Very few integrals, arising from diagrams with a charm 
quark in a closed fermion loop, depend in addition on the ratio z = mc/mb. We have 
suppressed this dependence in ([7]) for simplicity. 

Notice that an integral has different representations in terms of {S, V, n, m} because 
of the freedom to shift loop momenta in DR. It is the underlying topology, i.e. the inter- 
connection of propagators and external momenta, which uniquely defines the integral. 



{Qi,2) ^ {u{p')d{q,)u{q2)\QiMp)) 
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Figure 2: Full set of (naively) non-factorizable 2-loop diagrams. The bubble in the last four 
diagrams represents the 1-loop gluon self-energy. Only diagrams with final state interactions, 
i.e. with at least one gluon connecting the line to the right of the vertex with one of the upper 
lines, give rise to an imaginary part. 
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In the following we loosely use the word topology in order to classify the integrals. An 
integral with t different propagators Vi with rij > is called a t-topology. The integrals 
in the considered calculation have topology t < 6. 

The reduction algorithm makes use of various identities which relate integrals with 
different exponents {n,m}. The most important class of identities are the integration- 
by-parts identities [H] which follow from the fact that surface terms vanish in DR 



o omi cms 

d'kd'l — ^n^W^n. = 0, Ve{kJ}. (8) 

In order to obtain scalar identities we may contract ([8]) with any loop or external mo- 
mentum under the integral before performing the derivative. In our case of two loop and 
two external momenta we generate in this way eight identities from each integral. 

A second class of identities, called Lorentz-invariance identities [I2], exploits the fact 
that the integrals in ([7]) transform as scalars under a Lorentz-transformation of the ex- 
ternal momenta. In this way we may generate up to six identities from each integral 
depending on the number of external momenta. In our example with only two linearly 
independent external momenta p and q there is only one such identity given by 



p ■ q \ p^ q'^ -I- q^ p^ p^ g^' 

dp^^ dq^^ J dqi^ dp^' 



ami Qnis 

'^1 ---^^ = (9) 



In total we obtain nine identities from a given integral, each of the identities containing 
the integral itself, simpler integrals with smaller {n, m} and more complicated integrals 
with larger {n, m}. It is important that the number of identities grows faster than the 
number of unknown integrals for increasing {n,m}. Hence, for large enough {n,m} the 
system of equations becomes (apparently) overconstrained and can be used to express 
more complicated integrals in terms of simpler ones. Not all of the identities being linearly 
independent, some integrals turn out to be irreducible to which we refer as Mis. 

In the considered calculation we typically deal with systems of equations made of 
several thousands equations. The solution being straight-forward, the runtime of the 
reduction algorithm depends strongly on the order in which the equations are solved. As 
a guideline for an efficient implementation we have followed the algorithm from [13] . 

The reduction algorithm enables us to express the diagrams of Figure [2] as linear 
combinations of Mis which are multiplied by some Dirac structures. As the coefficients 
in these linear combinations are real, we may extract the imaginary part of a diagram at 
the level of the Mis which is a much simpler task than for the full diagrams. As depicted 
in Figure [3l we find 14 Mis that contribute to the calculation of the imaginary part of 
the NNLO vertex corrections. 



3.2 Calculation of Master Integrals 

Some Mis in Figure [3] can be solved easily e.g. with the help of Feynman parameters. For 
the more complicated Mis the method of differential equations in combination with 
the formalism of Harmonic Polylogarithms (HPLs) [TJ] turned out to be very useful. In 



6 



Figure 3: Scalar Master Integrals that appear in our calculation. We use dashed 
lines for massless propagators and double (wavy) lines for the ones with mass rrib 
(rUc)- Dashed/ solid/ double external lines correspond to virtualities / um'^ / m^, 
respectively. Dotted propagators are taken to be squared. 

this section we give brief reviews of these techniques and conclude with a comment on 
the calculation of the boundary conditions to the differential equations. 

The analytical results for the Mis from Figure [3] can be found in Appendix A.l of [Tj. 
As an independent check of our results we evaluated the Mis numerically using the 
method of sector decomposition ^16j . 

Method of Differential Equations 

The Mis are functions of the physical scales of the process which are given by scalar 
products of the external momenta and masses of the particles. In our calculation the Mis 
depend on the dimensionless variable m as in ([7]). 

For a given MI we perform the derivative with respect to u and interchange the order 
of integration and derivation 



The right-hand side being of the same type as equations ([H]) and (Q, this procedure 
again leads to a sum of various integrals with different exponents {n, m}. With the help 
of the reduction algorithm, these integrals can be expressed in terms of Mis which yields 
a differential equation of the form 




(10) 



d_ 

du 



Mli(M) = a{u]d)mi{u) + ^hj{u]d)Mlj{u), 



(11) 



where we indicated that the coefficients a and bj may depend on the dimension d. The 
inhomogeneity of the differential equation typically contains Mis of subtopologies which 
are supposed to be known in a bottom-up approach. In case of the Mis from the third 
line of Figure [3l one MI in the inhomogeneous part is of the same topology as the MI on 
the left hand side of ( ITTi) and thus unknown. Writing down the differential equation for 
this MI, we find that we are left with a coupled system of linear, ffist order differential 
equations. 

We are looking for a solution of the differential equation in terms of an expansion 

MI.H = J2'-^. (12) 

j 

Expanding (fTTj) then gives much simpler differential equations for the coefficients Cij which 
can be solved order by order in e. In addition, in the case where we were left with a coupled 
system of differential equations, the system turns out to decouple in the expansion. The 
solution of the homogeneous equations is in general straight-forward. The inhomogeneous 
equations can then be addressed with the method of the variation of the constant. This 
in turn leads to indefinite integrals over the inhomogeneities which typically contain 
products of rational functions with logarithms or related functions as dilogarithms. With 
the help of the formalism of HPLs these integrations simplify substantially. 



Harmonic Polylogarithms 

The formalism of HPLs [15j allows to rewrite the integrations mentioned at the end of the 
last section in terms of familiar transcendental functions which are defined by repeated 
integration over a set of basic functions. We briefly summarize their basic features here, 
focussing on the properties that are relevant for our calculation. 

The HPLs, denoted by H{fnw;x), are described by a w-dimensional vector rhw of 
parameters and by its argument x. We restrict our attention to the parameters and 1 
in the following. The basic definitions of the HPLs are for weight w = 1 

H{0;x)=\nx, H (1; x) = - \n{l - x) (13) 

and for weight w > 1 

H{a, m^-i, x) = dx' f{a; x') H{m^_i; x'), 
Jo 

where the basic functions /(a; x) are given by 



(14) 



/(O; x) = ^H{0; x) = -, /(I; x) = ^H{1; x) = (15) 

ox X ox 1 — X 



In the case of = 0, the definition in (1141) does not apply and the HPLs read 

H(0,...,0;x) = ^\irx. (16) 
wl 
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The HPLs form a closed and linearly independent set under integrations over the basic 
functions /(a; x) and fulfil an algebra such that a product of two HPLs of weight Wi and 
W2 gives a linear combination of HPLs of weight w = Wi + W2- 

As described above, the solution of the differential equations leads to integrals over 
products of some rational functions with some transcendental functions as logarithms or 
dilogarithms. More precisely, we find e.g. integrals of the type 

It turns out that all these integrals can be expressed as linear combinations of HPLs 
of weight w + 1. This is obvious for the first integral as it simply corresponds to the 
definition of a HPL, cf. (1141) with a = 1. For the other integrals in fll7p . an integration- 
by-parts leads either to a recurrence relation or again to integrals of the type f[T^ . Not 
all integrals in our calculation fall into the simple pattern f|T7|) . but a large part of this 
calculation can be performed along these lines. 

In the considered calculation we encounter HPLs of weight w < 3. Our results can 
be expressed in terms of the following minimal set of HPLs 

H{0;x)=\nx, H{0,0,1; x) = U^ix), 

H{1; x) = - ln(l - x), H{0, 1, 1; x) = Si,2(x). 

H{0,l;x) =U2{x), (18) 

The situation is more complicated for the last two Mis in the third line of Figure [3] where 
the internal charm quark introduces a new scale. However, a closer look reveals that 
these Mis depend on two physical scales only, namely uml and = z^ml. The Mis can 
then be solved within the formalism of HPLs in terms of the ratio ^ = /u if we allow 
for more complicated arguments of the HPLs as e.g. = | (l — a/1 + 4^) . 

Boundary conditions 

A unique solution of a differential equation requires the knowledge of its boundary con- 
ditions. In the considered calculation the boundary conditions typically represent single- 
scale integrals corresponding to m = or 1. It is of crucial importance that the integral 
has a smooth limit at the chosen point such that setting m = or m = 1 does not modify 
the divergence structure introduced in (fT2!) . 

In some cases the methods described so far can also be applied for the calculation 
of the boundary conditions since setting m = or 1 leads to simpler topologies which 
may turn out to be reducible. If so, the reduction algorithm can be used to express the 
integral in terms of known Mis. If not, a different strategy is mandatory. In this case 
we tried to calculate the integral with the help of Feynman parameters and managed in 
some cases to express the integral in terms of hypergeometric functions which we could 
expand in e with the help of the Mathematica package HypExp [T7]. Finally, the most 
difficult single-scale integrals could be calculated with Mellin-Barnes techniques [TS] . 
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4 Renormalization and IR subtractions 



The matrix elements {Qi) which we obtained from computing the 2-loop diagrams in 
Figure [2] are uhraviolet (UV) and infrared (IR) divergent. In this section we show how 
to extract the finite hard-scattering kernels T/ from these matrix elements. 



4.1 Renormalization 



The renormalization procedure involves standard QCD counterterms, which amount to 
the calculation of various 1-loop diagrams, as well as counterterms from the effective 
Hamiltonian. We write the renormalized matrix elements as 



(Qj)ren ^xl) ^ij {,Q j)ha.rei 



(19) 



where Z, 



ryl/2 ^3/2 



^ contains the wave-function renormalization factors of the 6-quark 



and Zq of the massless quarks, whereas Z is the operator renormalization matrix in the 
effective theory. We introduce the following notation for the perturbative expansions of 
these quantities 



(4) 



ren/bare 



oo , 



k=0 



ren/bare' 



oo 

E 

k=l 



An) 



k=l 



(k) 



(20) 



and rewrite f|T9l) in perturbation theory up to NNLO which yields 
/n.\(o) _ /n.\(o) 

Wj/ren W«/bare' 



(1) 
ren 



{Qi) 

{Qi)icn 



m 

(4) 



(1) , 

bare ' 

(2) , 
bare ' 



7(1) , 7(1)^ 



(4) 



(0) 
3 1 bare ' 



. (1) 
j I bare 



7(2) 7(1) A(l) 



7(2)^ 



(4-) 



i/bare- (^1) 



The full calculation thus requires the operator renormalization matrices For the 

calculation of the imaginary part, the terms proportional to the tree level matrix elements 
do not contribute and Z^^^ drops out in fl^ as expected for an effective 1-loop calculation. 

Mass and wave function renormalization are found to be higher order effects. For the 
renormalization of the coupling constant we use 



1 - 



as 
Aire 



11 

Y 



(22) 



according to the MS-scheme. The expression for the 1-loop renormalization matrix Z^^^ 
can be found e.g. in [19] and reads 



2 ^ 

^ 3 



_5_ 

12 



6 10 



1 

e' 



(23) 
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where the two hnes correspond to the basis of physical operators {Qi,Q2} and the four 
columns to the extended basis {Qi,Q2, Ei, E2} including the evanescent operators Ei^2 
defined in 



4.2 Factorization in NNLO 

In this section it will be convenient to introduce the following short-hand notation for 
the factorization formula ([1]) 

{Qi)ren = F ■ T, ® $ + . . . (24) 

where F denotes the B Mi form factor, Tj the hard-scattering kernels T/ and $ the 
product of the decay constant and the distribution amplitude 0m2- The convolution 
in ([1]) has been represented by the symbol (8> and the ellipsis contain the terms from 
spectator scattering which we disregard in the following. 
Formally, we may introduce the perturbative expansions 

^-E(S)^-". ^-E(£)^'". *-E(£)*'" (25) 

k=0 k=0 k=0 

Up to NNLO the expansion of (^M then yields 

mill = ■ rf ) + ■ ) ® $(0) + ■ Tf ) ® $(1), 

+ . ® $(0) + ■ Tf ) ® $(1) + ■ ® $(2). (26) 
In LO the comparison of ( 1211) and ( l26l) gives the trivial relation 

(4)(°)-(4)[,l = i^(°)-Tr®$(°) (27) 

which states that the LO kernels Tj^^^ can be computed from the tree level diagram in 
Figure H^. In order to address higher order terms we split the matrix elements into its 
(naively) factorizable (f) and non-factorizable (nf) contributions 

(4)£e - ^ + (28) 

The corresponding 1-loop diagrams are shown in Figure |Dd and Ht respectively. To this 
order (El and (ESI) lead to 



do) 
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(c) 







Figure 4: Tree level diagram (a), naively factorizable (h) and non-factorizable (c) 
NLO diagrams. 



which sphts into 



for the calculation of the NLO kernels t}^^ and 



(30) 



(31) 



which shows that the factorizable diagrams and the wave-function renormalization are 
absorbed by the form factor and wave function corrections F^^^ and 
This suggests in NNLO the following structure 



{Q^)?^ + z';'{Q^)r> + zl;>{Q,) 

= ■ t/°) ® + ■ T^'^ ® $(1) + ■ t/°) ® $(2). (32) 

These terms are thus irrelevant for the calculation of the NNLO kernels which justifies 
that we could restrict our attention to the non-factorizable 2-loop diagrams from Figure [2j 
In NNLO the remaining terms from 0211) and 0261) contain non-trivial IR subtractions 



^(1) 



(1) 



\(0) 



(1) 

nf 



(1) 



+ 



y{2) A(l) 



(33) 



This equation can be simplified further when we make the wave function renormalization 
factors in the form factor and the distribution amplitude explicit 



7I/2 yl/2 p 



z $ 



(34) 



Notice that the resulting amputated form factor -Famp and wave function $amp contain UV- 
divergences by construction. Using 03OI) . we see that the wave function renormalization 
factors cancel and arrive at the final formula for the calculation of the NNLO kernels T- 



(2) 



7(1) 



(1) 



\(0) 



= ■ Tf ) ® $(°) + F«p ■ T^'^ ® $(°) + ■ t/^) ® $«p- (35) 

As the tree level matrix elements and the factorizable 1-loop diagrams do not give rise 
to an imaginary part, these terms can be disregarded in the present calculation. 
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4.3 IR subtractions 



We now consider the IR subtractions on the right hand side of (135|) in some detaiL Let 
us first address the NLO kernels t/^"* which can be determined from equation ( l30l) . The 
renormahzation in the evanescent sector imphes that the left hand side of (1301) is free of 
contributions from evanescent operators up to the finite order e^. However, as the NLO 
kernels enter fl35|) in combination with the form factor correction F^mp which contains 
double (soft and collinear) IR divergences, the NLO kernels are required here up to 0{e^). 
Concerning the sub leading terms of 0{e), the evanescent operators do not drop out on 
the left hand side of (|30l) and we therefore have to extend the factorization formula on 
the right hand side to include these evanescent structures as well. Schematically, 



+ z 

nf ' • 



$(0) + i.{0) . T.d) 



.E 



(0) 



(36) 



with a kernel Tj^^ = (9(1) and an evanescent tree level matrix element F^^^^^^ = 0{ej^. 
Similarly, the right hand side of (!35|) has to be modified to include these evanescent 
structures 



(1) 



^(2), 



v(0) 



amp 



p(0) ^(2) ^(0) ^(1) , p 



amp 

(0) , p{0) ^(1) ^(1) 
E + ^E ■ J-i,E® ^amp,E- 



(37) 



(1') 

Notice that the term with the kernel ^ = 0{1) is not required to extract the finite 



(2) 

piece of the physical NNLO kernel . 

From the calculation of the 1-loop diagrams in Figure Ht, we find that the NLO 
kernels vanish in the colour- singlet case, Tg^'' = T^^^ = 0, whereas the imaginary part of 
the colour-octet kernels is given by 



1 



2lnu + 2lnu)(l + eL + =-e^L'^ 



+ 



;i 1 - 3 In M - In^ M + In^ m) (^e + e^L^ 
— - 26 + (^2 + yj Inu + (^9 - y ) Inw 



3 1 

- In^ M (in^ u - In^ u) 

2 3 ^ ' 



tin the notation of [2], the right hand side of (|36p corresponds to matrix elements of SCETj operators 
of the form [{^Wci)Tihy][{xWc2)T2iW^2X)]- In NNLO we match onto two SCETj operators with Dirac- 
structures Fi r2 given by Oi — rj.+L ® and O2 = r/l^^'^^'^L (g) ^jA_f^J_LijL (in our notation 
p' — ^mi,n_ and q = ^mftn^). The matrix element of Oi defines our structure _F('')<i)('^) and the 
evanescent combination 3O2 ~ 120i defines pj^'^^^jp . 
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1 



where L = In /i^/m^ and we recall that u = 1 



(38) 



u. 



Form factor subtractions 

We now address the form factor corrections which require the calculation of the diagram 
in Figure [5] (for on-shell quarks) and its counterterm. According to the definition of -Famp 
in 0341) . we do not have to consider the wave function renormalization of the quark fields. 

We again have to compute the corrections for physical and evanescent operators. 
Concerning the physical operator, the counterterm is found to vanish which refiects the 
conservation of the vector current. The evaluation of the 1-loop diagram in Figure [5] gives 



Fill $W = -C, 



amp 



2 

mi 



Tie] 



(39) 



which contains double IR singularities as mentioned at the beginning of this section. On 
the other hand, the 1-loop diagram with an insertion of the evanescent operator yields 
a contribution proportional to the evanescent and the physical operators. We now have 
to adjust the counterterm such that the renormalized (IR-finite) matrix element of the 
evanescent operator vanishes (which ensures that the evanescent structures disappear in 
the final factorization formula). We obtain 



^(1) ${o) 
-^amp.E 



mi 



rfe, Hl£(i±£. - 24 



ml 



(i-< 

1 -3£ 



i?(o) ${0) 



e(l -2e)(l-£)2 



(40) 



The form factor subtractions in ( !37l) then follow from combining (!38l) . ( |39l) and ( l40l) . 
We emphasize that the corrections related to the evanescent operator do not induce a 

(2) 

contribution to the physical NNLO kernel T| in this case since 



(1) ^(0) 



TT 



amp. 



0{e) 



$(o)_ 



(41) 



Wave function subtractions 

Concerning the wave function corrections we are left with the calculation of the diagrams 
in FigureOfor collinear and on-shell partons with momenta uq and uq. However, as in our 



q' = 



Figure 5: 1-loop contribution to the form factor correction F, 



(1) 

amp ■ 
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Figure 6: 1-loop contributions to the wave function correction <I>imp. The 
dashed line indicates the Wilson-line connecting the quark and antiquark fields. 

set-up = all these diagrams vanish due to scaleless integrals in DR. We conclude that 
the wave function corrections are determined entirely by the counterterms. We compute 
these counterterms by calculating the diagrams from Figure [6] with an (IR- finite) off-shell 
regularization prescription in order to isolate the UV-divergences (which are independent 
of the IR regulator). The counterterm for the physical operator is found to be 



2Cf 



dw V{u,w) F(°)$W(w) 



with the familiar Efremov-Radyushkin-Brodsky-Lepage (ERBL) kernel [20 

V(u,w) = \9(w-u)- ( 1 + ^^^ +9(u-w)- ( 1 + ^^ 
w \ w — u I w \ w — u 



(42) 



(43) 



where the plus-distribution is defined as [f{u,w)]_^ = f{u,w) — 6{u — w) Jq dv f{v,w). 
For the evanescent operator we obtain 



2Cf 



dw 



24:eVE{u,w) $(°)(w) + V{u,w) ^f{w 



(44) 



where Ve{u^w) denotes the spin-dependent part of the ERBL kernel given by 

11 u 

Ve{u,w) = 9{w -u)— + e{u-w) — . (45) 
w w 

Notice that the evanescent operators do induce a finite contribution to the physical NNLO 

(2) 

kernel in this case as the convolution with the corresponding NLO kernel implies 
1 



TT 



6C| 



+ Oie) 



i?(o)$(o)_ 



(46) 



We finally quote the result for the convolution with the physical NLO kernel 



-FWlmTf)<|.« = ^ 

_ 1 amp jY 



TT 



In M In M 



u 



u 



In^ M — 2 In M In M — In^ u — 4Lio iu 



1 ^\ TT 15 , 5m — 4 /In u Inn 

- + L +----2(3 + ^— — + 

e 2 A 2 \ u u 



21 

— mu 
3 



1 1 

— 3Li2(u) In^ u H 

^ ' 2u 2u 



3?/ i2— ^iS i2i — 

m u m u + m umu 

3 



+ ^ In^ u + 2 \nuU2{u) + 2Li3(M) + 2Si,2(?/) + 0(e) | F^") <1>(°). 

(47) 
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5 Vertex Corrections in NNLO 



We now have assembled all pieces required for the calculation of the NNLO kernels 
from (1371) . We indeed observe that all UV and IR divergences cancel in the hard-scattering 
kernels as predicted by the QCD Factorization framework. Since this is the result of a 
complicated subtraction procedure, this can also be seen as a very stringent cross-check 
of our calculation. 



5.1 ai in CMM basis 



The procedure outlined so far leads to the colour- allowed tree amplitude in the CMM 
operator basis defined in (jl]). We write 



ai(MiM2) =(72 + 



4^ 2Nr 



(4J 



where the ellipsis denote the terms from spectator scattering which are irrelevant for our 
purposes. In the CMM basis, the imaginary part of the vertex corrections V"(^'^) takes 
the form 



- Im f 

TT 



du gi{u) 0m2(m), 



- Im VI 



(2) 



du 



29 2 

—Ca - -nf)gi{u) + CFhi{u) 



In 4 

mt 



3 3 

+ Cp h2{u) + Ca h2.{u) + (n/ - 2) hi{u] 0) + /i4(-u; z) + hi{u] 1) \ 0a/2(-u), 



- Im K 

TT 



(2) 



du 



a 

Qgi{u) In — ^ + hQ{u 
mt 



(49) 



In writing (H9|) we have made the dependence on the renormalization scale explicit and 
disentangled contributions that belong to different colour structures. The NLO kernel 
gi{u) is given by 



-3 - 21nu + 21nM 



(50) 



and the NNLO kernels hi{u) will be specified below. The kernel hi{u\Zf) stems from 
diagrams with a closed fermion loop and depends on the mass of the internal quark 
through Zf = rrif/mf,. We keep a non-zero charm quark mass and write z = Zc = rric/mb 
for simplicity. 
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The NNLO kernels were so far unknown. They are found in this work to be 

ho{u) ■■ 



155 4 
^ + 4(3 + 4Li3(M) -4Si,2(m) - 12 In M Lis (m) + - In'^w - ein^wlnu 



H z — ^12(1*) — In H : Intilnti 



u 



-( 



u 



+ 



2u 
-^Inu — 



2uu 



4-nu + 2u^ , 47rS, (5 + 6x^2 - 12ii^)7r2 



241414 



+ (u ^ u) 



+ 



3-i4 + 7i/\ 2 ll-10i/^.. . . 1-Uu' 



+ 



2u 

46 - 51i4 

u 



In 14 — 
lni4 — 



4l414 

(41 - 42i4^)7r^ 
24^ 

2(13- 12m) 



U2{u) + 



4it 

— {u u) 



In ii In 14 



14 



hi{u) = 36 + 2 In^ 14 — 4Li2(i4) H — z In 14 — (14 <-> 14) 

g 

h2{u) = 79 + 32C3 - 16Li3(i4) - 32Si,2(^/) + - In^ 14 + 



2(4-1/2) 



u 



Li2(l4) 



13 - 9i/ + 6i/2 - 4m^ i7_6m2-8i/\ , 

in It H ^ in 14 in 14 



2i2 

'5-lli/ + 2i/2 47r2 



14 



+ 



In 14 



Auu 

(1 + Uu^ - Su^)tt^ 
8uu 



+ (14 <-> 14) 



+ 



2 9 — 1414^ 

4Li3 (it) + 4Si,2 (it) In^ It + 2 In^ It In 14 z — l^k{u) 

3 



uu 



13-1114+1414^ 2 

H In^ 14 + 

2ix 

24- 231/ 7i^\, 

z + — lni4- 

it 3 / 



5 - 714^ 



In It In It 



It 



+ 4 



(26 - 21u ^)n^ 
6u 



(it <-> 14) 



1 QTO 4 _ 7,2 

— - - I2C3 + 6U,{u) + 12Si,2(m) -In'u ^Li2(i4) 



+ 



+ 



4l4 



■ In It In It 



+ 



9 - 2i4 + 6142 - 4i4^ 2 7 + 414^-414^ 

m 14 z w " 

4l414 

, , (I + 6i/2 - 4i/4)7r2 ^ 
In 14 + z h (14 ^ 

' 8l414 

1. . 15-261/2 



/41 - 661/ + 8142 2\n (l + 6i/2-4i 
z h tt"' In 14 + -^^ z — 

\ 4l4 / 8l414 

-2Li3(i/) +4Si,2(m) +41ni/Li2(i/) + ^ln^i/ + 



+ (14 <-> 14) 

Li2(it) 



4lt14 



11 - 14i4 - 42i/2 , 2 11-141/2, , 
H 7^73 In It — In It In It 



12l4 



4l4 
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/2165 - 2156m 



TC 



y ) lllM + 



(53 - A2u^)tt^ 



24m 



(U ^ u 



hi{u] z) 



+ 



1! + !i _ 2e'ln2 ^ + ^lnz2 - (1 + 2^) Inn 

6 U X2 u 

+ (2{l + 4^X1 + 4^x2) Inxi - (Aixi + 2(1 + Ai)x^ \nx2 + {u ^ u) 

94^2 2(1 - ?>e 



9u 3 
12 + 29^ + 2^2 



:1-2m)(6m-<2) 
9uut, 



9e 



- m mulnz + 

X2 3 

liiM - ^(^(1 + 40(6 + 50x1 - 6(1 - 3e')x2 ) Inxi 



- ^(6(1 - 3e)xi - (1 + 40(6 + 50x2) \nx2 ~ (u ^ u) 



(51) 



where the last kernel has been given in terms of 



1 + 4^-1 
In the massless limit h4{u; z) becomes 

h^^u] 0) 



X2 = ^[^l+4i + l 



u 



17 2,2 2,2 - 20, 38, _ 

In u -\ — In u -\ mu mu. 

3 3 3 9 9 



(52) 



(53) 



5.2 ai and a2 in traditional basis 

Following our strategy from Section [2l we compute the colour-suppressed amplitude ^2 
by rewriting the colour-allowed amplitude ai in the traditional operator basis from ([3]). 
Manifest Fierz symmetry in this basis relates the two amplitudes via 



iWiiWo = H 



+ ... 

(54) 



where the upper (lower) signs apply for i = 1 (i = 2) and the elhpsis correspond to 
the terms from spectator scattering. In order to derive V^'^^'^^ we have to transform the 
Wilson coefficients in the CMM basis Ci into the ones of the traditional basis Cj. To 
NLL approximation this transformation can be found e.g. in ^ and reads 



a. 



14 ~ 



= 2^2 + ^ [AC^ + -C2 

C2 = C, + \C2 + ^{^C2 



(55) 
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Combining (HSj) . (13^ and we obtain 



- Im F(^) = - Im f 



TT 



TT 



gi{u) (pMiiu), 



- Im f/^) 

TT 



Im 



TT 



1 + 2 



Im 



(2) 



TT 



f du I -3gi{u)\n^ + 2gi{u) + l-ho{u)\ (f)M2{u 
Jo I ) 



Im 



/20 2 , 

(^yCA + QCf - -rif )gi{u) + Cphiiu) 



In At 

mt 



h2{u) - /io(m) + 45(1 (m) 



1 



/13(m) + -/Io(m) - gi{u) 



+ (n/ - 2) /;,4(u; 0) + /;.4(m; z) + /;.4(u; 1) > (pMii^). 



(56) 



Equation fISB]) represents the central result of our analysis, specifying the imaginary 
part of the colour-allowed tree amphtude ai and the colour-suppressed tree amplitude 
a2 according to (l5i|) . The expression for is in agreement with [1], whereas the 
expressions for VI2' are new. The kernels gi{u) and /^^(m) can be found in fl50l) and flSTj) . 
respectively. The terms proportional to n/ have already been considered in the analysis 
of the large /3o-limit in ^T]. Our results are in agreement with these findings. 



5.3 Convolutions in Gegenbauer expansion 

Our results in fl56|) have been given in terms of convolutions of hard-scattering kernels 
with the light-cone distribution amplitude of the emitted meson M2. We may explicitly 
calculate these convolution integrals by expanding the distribution amplitude into the 
eigenfunctions of the 1-loop evolution kernel 



(pA42 (u) = 6UU 



n=l 



(57) 



where a^^ and Ci^^^^ are the Gegenbauer moments and polynomials, respectively. We 
truncate the Gegenbauer expansion at n = 2 and perform the convolution integrals 
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analytically. Wc find 

/ du gi (u) (u) = -3-3a^\ 
Jo 

/ du ho{u) (j)M2{u) 
Jo 



1333 477r2 



15 237r2 



M2 



173 IStt^ 
+ 



30 35 

duhi{u) (pM^iu) = 36 + 28 af^ 



"2 ) 



/' 

^0 
^0 



du h2{u) (f)M2{u) = 



1369 1397r2 
~6~^^5 

103 7l7r2 



17 5l7r2 



+ 



du hz{u) (t)M2{u) 



15 35 

481 77r2 



,M2 
*2 ) 



,M2 



3 30 
1531 1697r2 



+ I2C3 - — + 



643 Utt^ 



12 10 



,M2 



80 



70 



,M2 

*2 ) 



{z)= I du hi{u] z) (t)M2{u) 
Jo 
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- + 148z' - 96^^F(z) - 36/ In' ^ 

3 ?/2 



l-(2yi + l)(l + 22z^) ln^-41ny2 



yi 



2/2 



A^ ,„2 2/1 



+ <^ 7 + 164z' + 180z^ + 144^" - 288^^F(z) + 12z*(3 + 16^^ + 12z^) In 



2/2 



1 - (2yi + 1)(1 + 22^2 + 84/ + 72/)] In^ - 41ny2 i of' 

J 2/2 J 

+ I - + 244/ + _ g4o^6 _ ggo^s ^ 24/(1 - 30/ - 40/) In' ^ 

15 3 y2 



576/F(^) + 8/(2|/i + 1)(6 + 11/ - 70/ - 120/) In ^ I af % (58) 



y2 



where we defined 



^ (^m^ - 1) , ?/2 = ^ (v^m^ + 1) , 



F{z) = Li3(-yi) - Si,2(-2/i) - lnyiLi2(-yi) + ^Inyiln't/s - ^In^-^' + (3- (59) 
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In the massless limit the function Hi{zf) simphfies to 

^4(0) = f + 7af^ + ^af . (60) 
6 5 

The finiteness of the convolution integrals in flSSl) completes the explicit factorization 
proof of the imaginary part of the NNLO vertex corrections. 

We summarize our results for the vertex corrections in the Gegenbauer representation 
of the light-cone distribution amplitude of the meson M2 (with Cf = ^,Ca = 3, Uf = 5) 

-Im \/W = -3-3af^ 

71 

i Im Vr = (9 + 9af") ,„ 4 + ™ + Hi! - 8C3 - (^ - 

71 ^ ^ ^ ^ ml 24: 90 ^ V 8 10 



,M2 



173 Qtt^ . 
60 ■ 35 ' ^ ' 



1 t7(2) _ fo. , HO.mA 10315 , 674vr2 28 



Im Vr' = - 26 + af ^ In ^ + Cs 

^ \ 3 ^ J ml 72 135 3 ^'^ 

10793 1667r2\ /3 1 55 1 877r2\ ^^^^ 



V 72 15 / ' V 48 42 



af^ + //4(^) + i^4(l), (61) 



with Hi{zf) given in (!58l) . In order to illustrate the relative importance of the individual 
contributions, we set = mh and 2; = mc/mb = 0.3 which yields 



Im V^(^) = -9.425 - 9.43af^ 

Im Vl^^ = 141.621 + 58.36 af^ - 17.03 af^ 

Im V^^^ = -317.940 - 115.62 af^ - 68.31 (62) 

We thus find large coefficients for the NNLO vertex corrections and expect only a minor 
impact of the higher Gegenbauer moments, in particular in the symmetric case with 
_ g_ Notice that all contributions add constructively in ai 2 due to the relative signs 
of the Wilson coefficients, C*i ~ 1.1 and C*2 ~ —0.2. In the case of ai the contribution 
from t^/^"* is found to exceed the formally leading contribution V^*^^-' due to the fact that the 
latter is multiplied by the small Wilson coefficient C2. Concerning 02 the NNLO vertex 
corrections are also substantial, roughly saying they amount to a 50% correction. In both 
cases the impact of the charm quark mass is small, we find a correction of ~ 3% compared 
to the massless case. A more detailed numerical analysis including the contributions from 
spectator scattering will be given in the following section. 

We finally remark that the large /3o-limit considered in [21] fails to reproduce the 

~ (2'] 

imaginary part of ai as it completely misses the leading contribution from Vi . In the 
case of a2 the approximation turns out to be reasonably good with a deviation of ~ 10% 
compared to the full NNLO result. 
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6 Numerical analysis 



6.1 Implementation of Spectator Scattering 

In the numerical analysis we combine our results with the NNLO corrections from 1- 
loop spectator scattering obtained in [21 [3], HI El |6] . In contrast to the vertex corrections 
considered in this work, the spectator term receives contributions from the hard scale 
fXh ~ "m-b and the hard-collinear scale ^hc ~ (AqcD^^fe)^''^- According to this, the kernels 
T/^ from ([T]) factorize into hard functions H-^ and a (real) hard-collinear jet-function J|[. 
Evaluating both kernels at the same scale fi would imply parametrically large logarithms 
which may spoil the convergence of the perturbative expansion. 

In order to resum these logarithms we follow reference [2] and perform the substitution 

Ciiflh) Hl\^h) ®l(ii{lJ,h,IJ,hc) ® Jwil^hc) ® [fB(pB]il^hc) ® (pMiifJ'hc) ® (pAhil^h), (63) 

where U\\ = U\\ consists of a universal Sudakov factor S and a non-local evolution 
kernel U\\. As an imaginary part is first generated at 0(«^) in the spectator term, 
we implement the resummation in the leading-logarithmic (LL) approximation. In the 
traditional operator basis from ([3|) the respective imaginary part takes the form 

Im a.iM,M,)\ = -)^^ QUfjEif^kc) y ^m,( ) ^m.( ) 



X Im 



+ 0(aD, (64) 



where we made the scale dependence of the parameters explicit and introduced the first 
inverse moment of the B meson light-cone distribution amplitude A^^. We further wrote 
1 = a^{n) in order to simplify the notation. In flMl) the resummation is encoded in 



y Jo Jo 



fj^h, fJ'hc) ri{u,z), (65) 



with the Sudakov factor S given in LL approximation in equation (106) of [6] and the 
kernels in equations (38) and (39) of [2]. Following [1] we defined 

Cmiz; fih, fJ'hc) = dv6vC^^'^\2v -l)U\\{v,z]fih,fJhc), (66) 



which can be computed by solving numerically the integro-differential equation 
d 



Cm{z; ij,, jj^hc) = - dw -f\\{z,w) Cm{w;fi,iJ.hc) (67) 
with initial condition Cmiz; fXhc, fJhc) = Cm {2z — 1) and 7|| from equation (99) of [6]. 
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In order to illustrate the numerical importance of the resummation we compare the 
values of the imaginary part of i?™" for m,n < 2 and fih = fJ'hc (line I, without resum- 
mation) and fih = 4.8GeV, fihc = l.SGeV (line II, with resummation) in Table [H We 
observe that the resummation leads to a suppression of the spectator term of ~ 10% due 
to the universal Sudakov factor (e""^ ~ 0.89 for our choice of input parameters). The 
resummation effects induced by U\\ turn out to be of minor numerical importance. 

According to flMl) we must evolve the Gegenbauer moments of the mesons Mi and 
M2 to the hard-coUinear and the hard scale, respectively. In LL approximation the 
Gegenbauer moments do not mix and the evolution reads 



Qs(/io) 



(68) 



64 and 72 = — 



with anomalous dimensions 71 — g ^^^^ — g 

We are left with the evolution of the B meson parameters to the hard-coUinear scale. 
We convert the HQET decay constant fs into the physical one using the LL relation 



/b(/^) 



a 



-2//3o 



(69) 



a^(mft) ) 

The evolution of is more complicated. The solution of the integro-differential equation, 
which governs the LL evolution of the B meson light-cone distribution amplitude, can be 
found in [22] . Here we adopt a model-description for the B meson distribution amplitude 
to generate the evolution of A^. We take the model from [22j which has the correct 
asymptotic behaviour and is almost form-invariant under the evolution. 

Finally we implement the BBNS model from |1] in order to estimate the size of power 
corrections to the factorization formula ([1]). This results in an additional contribution 
from spectator scattering related to subleading projections on the light-cone distribution 
amplitudes of the light mesons. It is given by 

_ -nasCp ^fMifs 

I power 



Im ai{MiM2)\ 



iV2 
c 



(70) 























. = 1 (I) 


11.0 


23.2 


29.4 


14.1 


23.8 


30.8 


15.1 


23.8 


31.3 


t = 1 (II) 


9.88 


20.9 


26.5 


12.5 


21.5 


27.8 


13.4 


21.5 


28.2 


^ = 2 (I) 


-5.29 


-8.43 


-8.24 


-6.58 


-11.0 


-11.3 


-7.04 


-12.0 


-12.6 


i = 2 (II) 


-4.77 


-7.60 


-7.43 


-5.85 


-9.72 


-9.98 


-6.27 


-10.6 


-11.2 



Table 1: Numerical values of the imaginary part of KP''^{ijLh-, l^hc) from fg5|] 
for fih = fific (line I, without resummation) and ji^ = A^.'&GeV, fi^c = l.5GeV 
(line II, with resummation). 
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Parameter 


Value 


Parameter 


Value 


A (5) 


u.zzo 


Jtt 


U.ioi 




4.8 


fR 
J B 


0.21 ± 0.02 


rric 


1.6 ±0.2 


F_f"(0) 


0.25 ±0.05 


rfibirhb) 


4.2 


AB(lGeV) 


0.48 ±0.12 


(m„ + md)(2GeV) 


0.008 ±0.002 


a^(lGeV) 


0.25 ±0.2 



Table 2: Theoretical input parameters (in units of GeV or dimensionless) . 



where r^{n) = 2m\j /mb{fi) / {rrig+mq) (fi) , Am = 1±^„(— l)"a*^ and X/^ parameterizes 
an endpoint-divergent convolution integral. The latter is written as 

XH = {l + PHe'^")ln'^ (71) 

Aft 

which may generate an imaginary part due to soft rescattering of the final state mesons. 
We take InniB/Ah — 2.3, pn = ^ and allow for an arbitrary phase (pn- 

6.2 Tree amplitudes in NNLO 

The numerical implementation of the vertex corrections from (15^ is easier since they 
can be evaluated at the hard scale ph ~ and depend only on few parameters. As can 
be read off from (l6Til . the first Gegenbauer moment of the meson M2 enters the leading 
term V'^^\ We therefore require its next-to- leading logarithmic (NLL) evolution (which 
can be found in [23]) but as we restrict our attention to i? — > tttt decays in the following 
the first moment does not contribute at all. Since the second Gegenbauer moment does 
not enter V'^^\ it is only required in LL approximation given by fl68|) . 

Our input parameters for the B ^ titt tree amplitudes are summarized in Table [2J 
The value for the B meson decay constant is supported by QCD sum rule calculations [21] 
and recent lattice results [25]. The form factor F^"^ (at large recoil) has been addressed 
in the light-cone sum rule (LCSR) approach [26j. As we implement the model from [22] 
for the B meson distribution amplitude, we take the respective value for A^. This value 
is somewhat larger than the one used in previous QCD Factorization analyses [H [21 El S] , 
but it is supported by a QCD sum rule and a LCSR calculation [27j. The value for the 
second Gegenbauer moment of the pion can be inferred from a LCSR analysis [28] and 
lattice results [29] . 

In order to estimate the size of higher-order perturbative corrections we vary the hard 
scale in the range ph = 4.8^24 GeV and the hard-collinear scale independently between 
Phc = l-S^oJ GeV. Throughout we use 2-loop running of with = 5 (n^ = 4) for 
quantities that are evaluated at the hard scale ph (hard-collinear scale phc)- The quark 
masses are interpreted as pole masses except for those entering . 
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With these input parameters the complete NNLO result for the imaginary part of the 
topological tree amplitudes is found to be 

Im ai(7r7r) = 0.012|^(i) + 0.03l|^(2) - 0.012|^(2) 

= 0.031 ± 0.015 (scale) ± 0.006 (param) ± 0.010 (power) 
= 0.031 ±0.019, 

Im a2(7r7r) = -0.077|^(i) - 0.052|^,2) + 0.020|^(2) 

= -0.109 ± 0.023 (scale) ± 0.010 (param) ± 0.045 (power) 

= -0.109 ±0.052. (72) 

In these expressions we disentangled the contributions from the NLO (1-loop) vertex 
corrections V^^\ NNLO (2-loop) vertex corrections V^'^^ and NNLO (1-loop) spectator 
scattering S^'^\ In both cases the NNLO corrections are found to be important, although 
small in absolute terms. For the imaginary part of ai the NNLO corrections exceed 
the formally leading NLO result which can be explained by the fact that the latter is 
multiplied by the small Wilson coefficient C2, cf. the discussion after ( l62i) . We further 
observe that the individual NNLO corrections come with opposite signs which leads to a 
partial cancellation in their sum. The phenomenologically most important consequence 
of our calculation may be the enhancement of the imaginary part of the colour-suppressed 
tree amplitude a2- 

In our error estimate in ( 1721) we distinguished between uncertainties which originate 
from the variation of the hard and the hard-collinear scale (scale), from the variation of the 
input parameters in Table [2] (param) and from the BBNS model which we used to estimate 
the size of power-corrections (power). By now we expect the power corrections to be the 
main limiting factor for an accurate determination of the amplitudes. However, although 
the inclusion of NNLO corrections has reduced the dependence on the renormalization 
scales, it still remains sizeable (in particular for ^h). For our final error estimate in the 
third line of each amplitude, we added all uncertainties in quadrature. 

Finally we remark that our numerical values for the NNLO spectator terms are much 
smaller than the ones quoted in [2J . This is partly related to the fact that we use different 
hadronic input parameters, in particular a much larger value for A^. In addition to this, 
the authors of [2] essentially evaluate the hard functions H-^ at the hard-collinear scale 
in order to partly implement the (unknown) NLL resummation of parametrically large 
logarithms. The NLL approximation is indeed required for the real part of the amplitudes, 
but as long as we concentrate on the imaginary part it is consistent to work in the LL 
approximation as discussed in Section 16.11 As the spectator term is the main source for 
the uncertainties from the hadronic input parameters, we also obtain smaller error bars 
than [2]. 
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7 Conclusion 



We computed the imaginary part of the 2-loop vertex corrections to the topological tree 
amplitudes in charmless hadronic B decays. Together with the 1-loop spectator scattering 
contributions considered in [21 [Sj HI El Ej, the imaginary part of the tree amplitudes is 
now completely determined at NNLO in QCD Factorization. 

Among the technical issues we showed that soft and collinear infrared divergences 
cancel in the hard-scattering kernels and that the resulting convolutions are finite, which 
demonstrates factorization at the 2-loop order. In our numerical analysis we found that 
the NNLO corrections are significant, in particular they enhance the strong phase of 
the colour-suppressed tree amplitude 02- Further improvements of the calculation still 
require a better understanding of power corrections to the factorization formula. 

Our calculation represents an important step towards a NNLO prediction of direct CP 
asymmetries in QCD Factorization. As the topological penguin amplitudes, which also 
affect the direct CP asymmetries, have not yet been computed completely at NNLO (the 
contribution from spectator scattering can be found in |4j), we refrained from discussing 
them already in this work. Moreover, the strategy outlined in this work may also be 
applied for the calculation of the real part of the topological tree amplitudes, which is, 
however, technically more involved f7\ IBU]. 
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